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ABSTRACT 

The model first introduced by Raga & Canto (1996) in which astrophysical jets are 
deflected on passing through an isothermal high density region is generalised by taking 
into account gravitational effects on the motion of the jet as it crosses the high density 
cloud. The problem is also generalised for relativistic jets in which gravitational effects 
induced by the cloud are neglected. Two further cases, classical and relativistic, are 
discussed for the cases in which the jet is deflected on passing through the interstellar 
gas of a galaxy in which a dark matter halo dominates the gravitational potential. The 
criteria for the stability of jets due to the formation of internal shocks are discussed. 
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1 INTRODUCTION 

According to the standard model for Fanaroff-Riley type II 
radio sources, a pair of relativistic jets of electrons is ejected 
in opposite directions from a v ery compact reg ion in the 
nucleus of an elliptical galaxy ( Blandford 199C ). The jets 
expand adiabatically through the interstellar medium of the 
host galaxy and the surrounding intergalactic medium. Col- 
limation of the jets is provided by the presence of a cavity or 
"cocoon" surrounding the jet, which maintains the pressure 
of the intergalactic medium in balance with that of the rela- 
tivistic gas within the jet. In this model, the jets maintain a 
straight trajectory. However, observations have shown that 
jets often bent. 

Dramatic changes in morphology are found in radio 
trails sources (for example NGC 1265), in which the jets 



beam passes through a region of interstellar or intergalactic 
gas in which there is a significant pressure gradient. This 
might occur if the jets pass through the interstellar medium 
of a large galaxy or diffuse intergalactic cloud. There are a 
few cases in which this may have occurred. In the case of 
3C34 Best, Longair & Rottgering (1997) found evidence that 



the jet had interacted strongly with the interstellar medium 
of a galaxy which happened to lie in the path of the jet. More 
evidence for misalignments between one of the hotspots of 
a radio galaxy and the host galaxy was found in 3C324 and 
3C368 ( Best et al. 1998 ). In these objects, no galaxy has been 
found in the path of the jets, but there is strong evidence 
for the presence of high density regions from the optical and 
infrared observations. For sources in which the jet has given 
rise to optically aligned emission regions, the optical emis- 



have a iemicircular shape, with the host galaxy at the pole. 
The curvature ot the jet is attributed to the motion ol the 
host galaxy through the intracluster medium which results 



sion is most probably produ ced by shock emission (Best 
Rottgering and Longair 200Cl[). The present study was mo- 



in a significant ram pre ssure acting on the jet (Begelman 
' Rees 198^ ). 



Blandford 



Some deflections may result from a combination of ki- 
netic and geometrical effects. For example, the precession of 
jets about a defined axis can cause the jet to be curved, al- 
though any fluid element in the jet always follows a straight 
trajectory. Another example is the proper motion of a galaxy 
through the intergalactic medium produc ed by the gravita- 



tional influence o f a companion galaxy (Begelman, Bland- 



ford & Rees 1984). Again, each fluid element in the jet fol- 



lows a straight trajectory. However, since the galaxy pro- 
ducing the jet moves in an orbit about its companion, the 
jet appears to be curved. 

Another way of inducing deflections in radio jets is if the 



tivated by these observations, the objective being to find a 
simple prescription for how significant the deflections of the 
radio jets could be if they passed through the interstellar 
medium of an intervening galaxy or gas cloud. In addition, 
we have analysed the stability of these deflections against to 
the formation of internal shocks. This stability analysis will 
be published elsewhere, but we discuss the results at the end 
of this article. 

The work presented here is an ext ension of the model 
first p roposed by Canto & Raga (1996) and Raga & Canto 
( 1996| ) in order to describe how defiections of radio jets are 
produced in jet/cloud interactions. The steady interaction 
of a jet with an isothermal sphere is analysed in two ways. 
Firstly we describe the steady interaction of a classical jet 
(i.e. a jet for which its internal velocities and thermal effects 
are not relativistic) with an isothermal cloud in which the 
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self gravity of the cloud is taken into acco unt. This is simi- 
lar to the analysis of Raga & Canto (1996), but generalised 



in order to take into account the effects of gravity. We also 
analyse the steady interaction of a relativistic jet with an 
isothermal sphere in which the gravity of the cloud is not 
taken into account. In this case, it possible to integrate the 
relevant equations exactly. Secondly, we consider a more rel- 
evant case for galaxies, in which it is assumed that all the 
gas in the galaxy is in hydrostatic equilibrium within a dark 
matter halo. Again, the cases of classical and relativistic jets 
are discussed. 



where e^, and are unit vectors in the directions x, 
r and respectively. Because angular momentum is con- 
served and the motion is two dimensional, the velocity is 
most simply written using eq.(^) as: 

V = VrSr -\ «osmv3oe,p, (3) 

r 

in which Vr — dr/dt represents the velocity in the radial 
direction. 

Since the steady flow of inside the jet expands adiabat- 
ically, we can calculate its path of by means of Bernoulli's 
equation: 



2 JET/CLOUD INTERACTIONS 

The interaction of an astrophysical jet with a cloud, in which 
the characteristic size of the cloud is much greater than the 
jet radius, has been studied in its initial stages by Raga 
& Canto (1995). They showed that initially the jet slowly 
begins to bore a passage into the cloud. Eventually a sta- 
tionary situation is achieved in which the jet penetrates the 
cloud and escapes from it in a direction which is different 
from the original jet trajectory. 

Once the steady state is reached, the trajectory is de- 
termined by maintaining pressure equilibrium with the sur- 
rounding environment. In other words, as the material in 
the jet moves, it adjusts its pressure in such a way that it is 
in equilibrium with the internal pressure of the cloud. Since 
the expansion of the jet is adiabatic, and a steady state has 
been reached, Bernoulli's equation can be used to describe 
the trajectory of the jet as it passes through the cloud. In 
what follows the interaction of a jet with a cloud of charac- 
teristic radius much greater than the diameter of an adia- 
batic jet is analysed, once the steady state configuration has 
been reached. 



3 CLASSICAL ANALYSIS 

In the classical case, the motion of a fluid element in the jet 
is determined by Euler's equation: 



dv 

— — h V ■ grad V 

at 



—grad p — grad ( 

P 



(1) 



where v is the velocity of the flow, and p and p are pres- 
sure and density. The gravitational force produced by the 
cloud is described in the term which contains the gravita- 
tional potential (j>. The left hand side of eq.(|l|) is the force 
per unit mass experienced by a fluid element as it moves. 
Since all quantities depend only on the distance from the 
cloud's centre r, vector multiplication of the radius vector r 
with eq.(^) implies rxdv /dt = Q. In other words, the speciflc 
angular momentum l = rxv is conserved as the fluid moves. 
Since the radius vector is perpendicular to the angular mo- 
mentum vector, the motion is two dimensional, and so, polar 
coordinates {r,ifi) are used in the following analysis. 

Consider a situation in which the jet enters the cloud 
parallel to the x axis at a distance ro, so that its velocity 
vector is initially given by: 



Vo 



'VOGx = —VO (cos ifo Br — siu ipo B^) , 



(2) 



i y dw^ + j dw + j 



d(/> = 0. 



(4) 



in which the line integrals are taken from the initial position 
of a given fluid element to its final one. The enthalpy per 
unit mass w of the flow in the jet is given hy w — T~^p/ p 
and r= (7 — l)/7 for a gas with polytropic index 7. 

Substituting eq.(|^ and eq.(Q) into the first integral of 
eq.(0) gives: 



dv 



ro vq sin ipo dr 
r r dip 



1 2 
+ Vo 



ro_y 
r / 



sin (/3o — 1 



(5) 

where we have used the fact that along the jet trajectory 
dr /vr = rd(p /v^. Since the gas in the jet obeys a polytropic 
equation of state, we obtain for the second integral in eq. 



7r Vpo 



(6) 



in which Cq is the initial sound speed of the jet material. 
The integral for the gravitational potential produced by the 
self-gravitating cloud is obtained from: 



d<t)^G 



M{r) , , ^ f dr 
dr = A-kG / -5- 



rpc(Ode, 



(7) 



where M(r) and pc(r) represent the mass and density of the 
cloud at a distance r. Substitution of eqs.(|5|)-(^ into eq.(^ 
gives the equation for the path of the jet as it expands: 



dp 



sin po 



rj^ sin^o 



7rM2 



- 1 



SttG f dr 



1/2 



(8) 



in which r] = ro/r and Mo is the initial Mach number of 
the flow in the jet. The positive and negative choices for the 
value of the derivative drj/dp in eq.(H) have to be chosen 
with care. For example, if we consider the case in which no 
gravity and no pressure gradients are taken into account (i.e. 
the last two terms on the right hand side of eq.(^) are zero, 
which corresponds to a straight trajectory) the derivative 
dr;/d<^ < for r;, < l/siniySo and vice versa. The equality 
r;, = l/sin(/Po corresponds to the point at which a given 
fluid element in the jet reaches the y axis during its motion 
in this particular case. 
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In the limit of high initial supersonic motion (Mo 3> 1) 
the third and fourth terms on the right hand side of eq.(p|) 
are important only when 77 = 1/ sin (po and we can simplify 
eq.(^ by making an expansion about this point. In general 
terms, if: 



— ] = a + /3?7 sin (^0 + (v^ sin Vo , 



and 



— 47rG [a + f3ri sin 1^20 + (rj^ sin^ifo 



(9) 



(10) 



are expansions of the pressure and gravitational potential 
respectively around ri=l/ sin(^o, then eq.(0) takes the form: 



dr? ^ _^ 
dip 

in which: 



sin ifo 



{a + br] + erj } 



211/2 



(11) 



, 2 (a - 1) 8ttG _ 



( 2/3 



+ 



SttG - 



13 sinipo, 



e=-\l + 



2C 



+ 



SttG 



C sin ipo, 



to second order in rj sin ipo. For the cases considered in this 
analysis, the general solution of eq.(|llj) is: 



?7 = — < V A sin 



+ 



yo — y 
sin ipo 

2e + b 



(12) 



with A = b^ — 4ae. The angle between the velocity vector of 
the jet streamline and the x coordinate axis on its way out 
of the cloud (deflection angle) can be calculated from the 
relation taxitp = {vy / v.x)cxit, in other words. 



ip = arctan 



sin (fc {drj/dif)^ — cos (pc 



+ 7r. 



_ cos (fie (dri/ dip)^ + sin pc 
where the exit azimuthal angle pc is given by: 



sm I „ . 2e + b 
ipe — po -\ ; •{ 2 arcsm — — h tv 



(13) 



(14) 



provided the deflections are not large. The derivative 
{d'ri/d'p>)c is evaluated at 77 = 1, with a negative choice of 
sign in eq.(^. 

3.1 Isothermal Cloud 

Let us consider the case of an isothermal cloud, for which the 
density in the cloud pdoud = (./r^ ■ In other words, because 
the jet and cloud are maintained in pressure balance, the 
pressure on the jet is given by: 



P_ 

Po 



To 



(15) 



For this isothermal case, it is easy to verify that: 



13-- 



1 - F (3 - 2r) 

sm (po 

4r (1 F) 

9r ' 
sm ipo 

F (2F - 1) 
sm (fio 



a = ,tln(sin(^o) + 



/3 = 2e, 



(16) 



c 



(17) 

Th is sol ution corresponds to that found by Raga & 
Canto (1996) for the case in which no gravity is present, 
i.e. a = f3 — ( — 0. From the solutions obtained above in 
eq.(p4|) and eq. (p^) it follows that the dimensionless quan- 
tity A = G^/MoCo is a number that parametrises the re- 
quired solution. 

The deflection of jets in isothermal clouds may be im- 
portant for interstellar molecular clouds and the jets associ- 
ated with Herbig-Haro objects. For this case we can obtain 
a value for the parameter A. If we adopt a particle number 
density of nn ~ 10^ cm^'^, and a temperature r~ 10 K for a 
cloud with radius ro ' 
then 



Ipc (3pitzerl99S; Hartmann 1998) 



A 



10" 



1 pc 



no 



T 



10" 



T 



100 kpc (Longair 1992; Longair 1998). With 



102cm-3; VlOK/ ■ ^"^^^ 

The same calculation can be made for the cases of radio 
jets interacting with the gas inside a cluster of galaxies. 
For this case, typical v alues are nu 
10^ K and ro ' 

these values, the parameter A ~ 10~^/Mo, exactly as 
eq.(^). The parameter A is an important number which 
can be obtained by dimensional analysis. For, the problem 
in question is characterised by the gravitational constant 
G, a "characteristic length" ro and the values of the veloc- 
ity of the jet and the density vo and po respectively. Three 
independent dimensions (length, time and mass) describe 
the whole hydrodynamical problem. Since four independent 
physical quantities (G, po, vo, andro) are fundamental for 
the problem we are interested, the Buckingham II-Theorem 
of dimensional analysis requires the existence of only one 
dimensionless parameter A = G^/M§Co, exactly as in the 
analytic solution. The fact that jets are formed in various 
environments such as giant molecular clouds and the gaseous 
haloes of clusters of galaxies with the same values of the di- 
mensionless parameter A provides a clue as to why the jets 
look the same in such widely different environments. 

From its deflnition, the parameter A can be rewritten 
as A = (3/47r)(GM/r)(l/iio), where M is the mass within 
a sphere of radius ro. This quantity is roughly the ratio of 
the gravitational potential energy from the cloud acting on 
a fluid element of the jet, to its kinetic energy at the initial 
position ro. The parameter A is thus an indicator of how 
large the deflections due to gravity are going to affect the 
trajectory of the jet. The bigger the number A, the more 
important the deflection caused by gravity will be. In other 
words, when the parameter A 1 the jet becomes ballistic 
and bends towards the centre of the cloud. When A ^ 1 the 
deflections are dominated by the pressure gradients in the 
cloud and the jets bend away from the centre of the cloud. 

Fig.(^ shows plots for three different values of A with 
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initial Mach numbers of Mo = 5 and A/o = 10. A compari- 
son with a numerical integration of eq.(|Tl]) using a fourth- 
order Runge-Kuta method is also presented in the figures by 
dashed lines. This comparison shows that as long as the de- 
flections are sufficiently small, or as long as the Mach num- 
ber of the flow in the jet is sufficiently large, the analytic 
approximations discussed above are a good approximation 
to the exact solution. 



3.2 Gas on Dark Matter Halo 

Let us consider next the case of a galaxy dominated by a 
dark mat ter halo for which the ma ss density is given by the 
relation (Binney & Tremaine 1987): 



The dimensionless number k which parametrises the 
required solution can be obtained by dimensional analysis 
(apart from an unimportant factor of — 47r) in the same 
way as the parameter A was obtained in Section 3.1. In 



this case the important parameters in the problem are the 
gravitational constant G, the characteristic length a and 
the sound speed c*. The number k can be rewritten as 
k = -(4/(4-7r))(GM/a)(l/cJ), where M is the mass of 
a sphere with radius a. This quantity is proportional to the 
gravitational energy of the cloud evaluated at the core ra- 
dius divided by the sonic kinetic energy that a fluid element 
in t he jet has. In other words, in the same way as in Section 
3.1, the dimensionless number k is an indicator as to how 



big deflections produced by gravity are. 



Pd = 



1 + 



(19) 



in which a is the core radius and quantities with a star refer 
to the value at the centre of the galaxy. 

The potential resulting from such a density profile can 
be calculated by means of eq. (^) , 



4 RELATIVISTIC ANALYSIS 

Let us consider the case in which relativistic effects are in- 
cluded in the interaction of a relativistic jet and a stratified 
high density region. In order to simplify the problem, the 
self gravity of the cloud acting on the jet is ignored. For this 
case, the relativistic 4-Euler's equation is: 



0d - = 4-KGpd^a X 



■In 



H — arctan 
r 



(20) 



in which the value of the gravitational potential (j}d^ has been 
evaluated at the centre of the galaxy r* = 0. If the gas in the 
galaxy is in hydrostatic equilibrium with the dark matter 
halo, then gradp= — pgrad 0^. In this case, the enthalpy 
of the isothermal gas is given by: 



t>d + (pd, = In ( — 



and so the pressure is: 

p = p^ exp{(-(^d + 0dj/co} 



(21) 



(22) 



It is possible to simplify the above expressions using the 
fact that for astronomical cases ro':S>a- In other words: 



P -k 
— = V , 
Po 



(23) 



where k = —ATvGpd^a^/ci- In this approximation, the re- 
quired analytic solutions can be found: 



a = i (rfc + 1) (Lfc + 2) sin^Vo, 
/3 = -rfc(2 + rfc)sin^'=<^o, 
C = irfc(rfc + l)sin^'=<y5o. 



(24) 



where for simplicity it was assumed that a = {3 = C^=0, in 
other words we neglect the gravitational field induced by 



the m ass of the cloud. Typical values (Binney & Tremaine 



1987) for galaxies are ~O.lM0pc , a ~ 1 kpc. Taking 



central values for the gas in the galaxy as ~ 1 cm and 
T*~10^K then fc~-10. Fig.(|) presents plots for different 
values of Mo and k. 



(25) 



k du, dp k dp 

ox" ox^ ox" 
in which w, and p are the enthalpy per unit volume and 
the pressure of a given f luid element in its prop er frame 
of reference respectively (Landau & Lifshitz 1987). The 4- 



velocity u ' = dx /ds (latin indices run as 0, 1, 2, 3) is given 
by (7,7w/c), c is the speed of light and 7 is the Lorenz 
factor. The 4-radius vector x'^ = [ct, r) and the metric tensor 
rjki = diag(l, — 1, — 1, — 1). The space component of eq.(|25|) 
can be written as: 



w"/~ I dv , 
' — — h V ■ grad V 

ot 



V dp 

-gradp^--. 



(26) 



which is the generalisation of Euler's equation. As before, the 
term in brackets on the left-hand side represents the classical 
force per unit mass acting on an element of fluid as it moves. 
By considering steady flow and noting that all quantities 
depend only on the radial coordinate r, vector multiplication 
of the radius vector r with eq.(p^ shows that the quantity 
I — r X V is conserved as the fluid moves. This quantity 
corresponds to the classical equivalent of the speciflc angular 
momentum, but it is not the relativistic angular momentum 
which is given by r x 71;. The constancy of I shows that the 
motion is two dimensional and so eqs.(g)-(^) are valid in the 
relativistic case as well. 

The trajectory of the jet is described by Bernoulli's 
equation, which is given by the time component of eq.(|25[) 
for steady adiabatic flow by (Landau & Lifshitz 1987): 



0. 



(27) 



The line integral is taken from the initial position of a 
given fluid element to its final one. The particle number per 
unit proper volume is n and we assume that the electrons 
in the jet are ultrarelativistic, so that the equation of state 
is given by p = e/3 with e being the internal energy density 
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Figure 1. Deflection produced in a jet due to the interaction with an isothermal cloud (semicircle) of radius tq. The jet penetrates the 
cloud from the right, parallel to the x axis. Different trajectories are shown in each diagram for different initial heights of y/ro =0.05, 
0.15,. ..,0.95 as measured from th e x/rg axis. In each figure the top diagram corresponds to the case in which gravitational effects are not 
considered (Raga & Canto 1996). The central and bottom diagrams represent trajectories for which gravitational effects are taken into 
account and the parameter A has values of 10~®, and 0.01 respectively in units of the square of the initial Mach number Mg of the jet 
(see text). A polytropic index 7 = 5/3 for the flow in the jet was assumed. The diagrams at the left and right correspond to initial Mach 
numbers for the jet flow of Mq = 5 and Mg = 10 respectively. The dashed lines in the graphs represent the direct numerical integration 
of the equation of motion. The continuous lines are the analytic approximations discussed in this article. 



of the fluid. The requirement that the pressure of the jet 
equals that of the cloud, together with the fact that pccn'^^^ 
gives: 



d?7 ^ ^ c 

dip vo sin (po 



, '"o 2 ■ 2 
1 -T] sm (po 



-2 I P 



1/2 



1/2 



(28) 



in which ri = ro/r. The sign of dri/d<p in eq.(|28[) varies as the 
jet crosses the cloud. For example, for the ultrarelativistic 
case, in which a straight trajectory is expected, it is pos- 
itive for 77, > 1/ sin ipo , and negative when the inequality 
is inverted. A general analytic solution of eq.(|2^) can be 
found by noting that for high relativistic velocities the third 
term on the right hand side of eq.(^8|) is important only for 
r] = 1/ sini^o- In other words, the pressure stratification of 
the cloud can be written as eq. (|) with r 1/2. We can 
therefore expand eq.(p^) around 77 = 1/sin^o to obtain a 
relation which is the same as eq.([Lll) but with: 



, b- 

7o 



f3 



2 smipo, e - 



+ 1 ) sinVo. (29) 



to find an exact solution to the problem, since eq.(]l5|) and 
eq.di^) imply: 



d?? , 1 / -2n 2 . 2 1 

:r- = i — (1 - ^770 ) ~n sm <p>o I 

dip sm (po 



(30) 



In other words, the solution is the same as the one already 
found in eq.(pT[) and eq.(p^ but with: 



«o 



^ ^ 7 ^, e = - sin^ Po. (31) 
Wo/ 



Fig.(|3|) shows plots of the trajectory of the jet for different 
values of the initial velocity of the jet. 



4.2 Gas on a Dark Matter Halo 

As in section for ro 3> a, the variation of the pressure 
in the galaxy is given by eq.(^) and the trajectory of the 
path of the jet is given by eq.(ll) and eq.(h2|) together with 
eq.(|24|) and the substitution F— > 1/2. Fig.(^ shows plots of 
this for fc = — 3 and difi'erent values of the initial velocity of 
the jet vq. 



4.1 Isothermal Cloud 



As in section 3.1, consider the case in which the jet propa- 
gates through an isothermal cloud In this case it is possible 



5 DISCUSSION 

When a deflection is observed, it is possible to work back- 
wards and flnd useful properties concerning the initial inter- 
action of a jet with a stratified density region. For example, 
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Figure 2. Deflection produced in a jet as it travels across a galaxy, for which its gravitational potential is dominated by a dark matter 
halo. The jet penetrates the galaxy parallel to the x/tq axis. Various trajectories are shown in each diagram for different initial heights 
= 0.5, 0.15,...,0.95 measured from the x/tq axis. The left and right diagrams were calculated for the case of a non relativistic jet with 
an initial Mach number of Mo = 10 and Mo = 20 respectively. For each of this diagrams a value of fc = — 1, —2, —3 was used for the too , 
middle and bottom panels respectively (see text). The dashed lines in the figures represent the direct numerical integration of eq.(|f l[) 
with the pressure given by eq.(E2|) for the case in which the ratio of the core radius a to the initial radius ro is given by a/ro = 10^'' <C 1. 
The continuous lines are analytic approximations found under this conditions. 



by taking the "standard" values used in sections 3.1 and 3.2 



for the pressure and density in the stratified gas it is pos- 
sible to calculate the initial azimuthal angle yn for a given 
initial velocity of the jet. Indeed, from eq.(|l3|) and eq.(^), 
together with the negative derivative (djy/dt/p)^ at the mo- 
ment the jet leaves the cloud we can find the value of the 
deflection angle cos t/j as a function of the velocity of the jet 
vq and the initial azimuthal angle <po- The curves for which 
cos ij) = const give the required relation between the initial 
velocity and azimuthal angle. Fig.(P) shows two examples of 
this type of analysis. 

Different combinations of the various parameters in- 
volved (or the known observables) in the problem can be 
formed so that, for a given deflection, the other quantities 
can be calculated. For instance one can ask for the values 
of the central density of the gas in the cloud, the density in 
the jet, etc. 

The assumption of pressure equilibrium between the 
jet and the external envi ronment has certain limitations on 
the geometry of the jet (Raga & Canto 1996). In the non- 



relativistic limit pressure equilibrium is achieved, when the 
time for a pressure scale height (p(r)/|dp(r)/dr|) to change 
signiflcantly within the jet material is greater than the sound 
crossing time of the jet radius; 



p(r) 



Wj \Ap{r)/dr\ 



(32) 



where cj is the sp eed of sound in the jet, nj its velocity 
and Aj its radius (Raga fc Canto 1996). For the case of an 



isothermal sphere the right hand side of eq.(|32[) is r/2vj and 
for the case of gas in pressure equilibrium with a dark matter 
halo is r/|fc|i;j. In the case of relativistic flow, eq.(B2|) is still 



valid but with the substitution v\ 



7„.«j and cj 7cjCj, 



where 7„- and 7cj are the Lorentz factors of the fiow velocity 
and the velocity of sound respectively. 

When a jet bends it is in direct contact with its sur- 
roundings and entrainment from the extern al gas might 
cause disruption to its structure (Icke 1991). However, if 
this situation is bypassed for example, by an efficient cool- 
ing, then there remains a high Mach number coUimated flow 
inside a curved jet. When supersonic flow bends, the char- 
acterist ics emanating from it int ersect at a certain point in 
space ( Landau fc Lifshitz 1987 ). Since the hydrodynami- 
cal values of the flow in a characteristic line have constant 
values, the intersection causes the different values of these 
quantities to be multivalued. This situation can not occur 
in nature and a shock wave is formed. 

The formation of internal shocks inside the jet gives rise 
to subsonic flow inside the jet and collimation may no longer 
be achieved. If the characteristic lines produced by the flow 
inside the jet intersect outside the jet, then a shock wave 
is not formed and the jet remains coUimated as it bends. 
However, two important points have to be considered in the 
discussion. The first is that the Mach number decreases in 
a bend as the fiow moves. The second is that the rate of 
change of the Mach angle with respect to the angle the jet 
makes with its original straight trajectory (the bending an- 
gle) increases without bound as the velocity of the fiow tends 
to that of the local velocity of sound. This was first proved 
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Figure 3. Deflection of a relativistic jet produced by its interac- 
tion with an isothermal cloud (semicircle). The jet is assumed to 
travel parallel to the x axis at the moment it enters the cloud from 
the right. In each plot different trajectories are shown for different 
values of the initial height of the jet y/rg =0.05, 0.15,. ..,0.95 as 
measured from the x/ro axis. The top, central and bottom panel 
plots were calculated for values of the initial velocity of the jet dq 
in units of the speed of light c of 0.99, 0.97, 0.95 respectively. 



by Icke (1991) for the case in which no relativistic effects 
were taken into account. We have made a relativistic gen- 
eralisation to these two points which will be discussed in a 
future paper. The increase in the Mach angle means that 
there has to be a stage at which this angle grows faster than 
the bending angle. This causes the characteristics in the jet 
to intersect and to form a shock. The only way to avoid this 
situation is if the rate of change of the Mach angle with re- 
spect to the bending angle is less than one. The formation of 
this shock is most probably not too destructive for the jet. 
This is because the Mach number of the flow is close one, 
which implies that the shock is weak. 

The fact that two shocks might form inside a jet, one 
of them at the end of the bending when the Mach number 
is near one, enab les us to find an upper limit to the bending 
angle (Icke 1991). For example, a classical jet with 7 = 5/3 



cannot be deflected more than 60.8°. Under the same classi- 
cal conditions, but by assuming a polytropic index 7 = 4/3, 
classical jets with a relativist ic equatio n of state cannot be 
deflected more than 148.12° (llcke 199l|). 



When a full relativistic analysis is introduced in this 
description, this upper limit can not be greater t han its 
classical counterpart. The relativistic Mach angle (Konigl 



198C) cecreases as the velocity of the flow approaches that 



of light. This means that characteristic lines emanating from 
a relativistic flow are closer to the streamlines as compared 
with their classical counterparts. The fact that characteristic 
lines are closer to the streamlines means that when a bend- 
ing occurs, the chance for an intersection between different 



Figure 4. Trajectory of a jet as it crosses a galaxy. The gas in 
the galaxy is assumed to be in hydrostatic equilibrium with a 
gravitational potential given by a dark matter halo in the galaxy. 
It is assumed that the jet enters the galaxy parallel to the x axis 
at a height of y/ro = 0.05, 0.15,...,0.95 in different cases. The plots 
were calculated for the case in which the parameter fc = — 3 and 
the initial velocity of the jet in units of the speed of light is 0.999, 
0.995 and 0.99 from top to bottom. Continuous lines are analytic 
approximations for a high Mach number of the flow inside the jet. 
The dashed lines are direct numerical solutions. 



characteristics increases. As a result, a relativistic jet with 
a polytropic index 7 = 4/3 can not bend more than 47.94°. 
The precise conditions under which an internal shock is pro- 
duced for a given jet depend on the shape of the curve that 
the jet makes as it bends and the radius of the jet. The 



e jets a 



are valid for 



general conclusion is that, provided the jets are sufficiently 
narrow, the deflections shown in 
small deflections, say < 20-30°. 

As an example, we c an take the case of the western je t 
in the radio galaxy 3C34 (Best, Longair & Rottgering 1997). 



If this jet was bent because of the interaction with a typical 
galaxy for which its gas is in pressure equilibrium with a dark 
matter halo, then we can use the standard values presented 
in section 4.2. First of all, the western jet in 3C34 has a 
bending angle Q ~ 10°, so that the deflection angle ~ 
170°. If we assume that the flow inside the jet moves with a 
velocity v = 0.99c, it follows from the right diagram in fig.(^ 
that sin Lpo ~ 0.3, or ipo ~ 17° . On the other hand, the value 
O ~ 10° is well below the upper limit of 47.94° discussed 
above, so that at least no terminal shock will be produced 
by the deflection of the jet. From preliminary results that 
will be published elsewhere, it follows that if the trajectory 
of the jet in 3C34 is circular, then in order not to produce an 
internal shock at the onset of the curvature, the ratio D/R 
{D is the width of the jet and R its radius of curvature) has 
to be less 0.08. In other words, the destruction of jets only 
occurs if the jets are thick, rather than narrow. This result 
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Figure 5. Variations of initial azimuthal angle i/pq as a function of velocity vq in units of the speed of light c for constant values of the 
deflection angle j/i. The angle ij) is the azimuthal angle the velocity vector of the flow in the jet makes with the x axis at the moment 
it leaves the cloud. Every plot was calculated for cosip = const with values given by i/) of 175° 170°...,155° . The gradient of ip decreases 
towards the lower left part of each diagram. In other words, deflections become stronger as the curves approach this region on the 
diagram. An isothermal sphere and an isothermal gas in hydrostatic equilibrium with a dark matter halo were assumed for the left and 
right diagrams respectively in the case of a relativistic jet. 



is in accord with the narrowing of jets observed in a wide 
range of different sources. 

The most important consequence of the calculations 
presented in this article is the sensitivity of the deflection 
angles to variations in velocity (see for example fig.(^)). This 
sensitivity is due to the fact that the force applied to a given 
fluid element in the jet (due to pressure and gravitational 
potential gradients) is the same independent of the velocity 
of the flow in the jet. However, as the velocity of the flow in 
the jet increases, there is not enough time for this force to 
change the curvature of the jet soon enough, giving rise to 
very straight jets. 
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